Potts model on complex networks 
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We consider the general p-state Potts model on random networks with a given degree distribution 
(random Bethe lattices). We find the effect of the suppression of a first order phase transition in 
this model when the degree distribution of the network is fat-tailed, that is, in more precise terms, 
when the second moment of the distribution diverges. In this situation the transition is continuous 
and of infinite order, and size effect is anomalously strong. In particular, in the case of p — 1, we 
arrive at the exact solution, which coincides with the known solution of the percolation problem on 
these networks. 
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I. INTRODUCTION 

Complex networks display a spectrum of unique effects 
[1-9]. Cooperative phenomena in complex networks are 
attracting much attention these days. The theoretical 
study of various cooperative models on random networks 
[10-23] has demonstrated that their critical behaviour is 
extremely far from that on regular lattices and on 'planar 
graphs' [24]. 

In this paper, we report our exact and asymptotically 
exact results for the thermodynamic properties of the p- 
state Potts model on uncorrelated random networks with 
a given degree distribution. These networks are the undi- 
rected graphs, maximally random (i.e., with maximum 
entropy) under the constraint that their degree distribu- 
tion is a given one, P(k). Here, degree is the number 
of connections of a vertex. Correlations between degrees 
of vertices in such graphs are absent, as well as cluster- 
ing. In graph theory, these networks are called 'labeled 
random graphs with a given degree sequence' or 'the con- 
figuration model' [25]. One should stress that this is the 
minimal model of complex networks. Most of results on 
cooperative models on networks were obtained just for 
this basic construction [18-21,23,26-29]. See, however, 
Refs. [30-33], where the Berezinskii-Kosterlitz-Thouless 
percolation phase transition was studied in growing net- 
works. See also Ref. [34] and references therein for co- 
operative models on networks with correlations between 
degrees of the nearest-neighbour vertices. 

The percolation problem on uncorrelated networks 
with a given degree sequence has been studied in 
Refs. [18-20]. The Ising model on these networks 
was studied by simulations in Ref. [10] and solved in 
Refs. [11,12] (see also Ref. [17]). It was shown that 
the presence of fat tails in the degree distributions of 
networks dramatically changes the critical behaviour of 
these models. But the site percolation problem and the 
Ising model are only particular cases of the general p- 
state Potts model [35,36]. The site percolation is equiv- 
alent to the one-state Potts model [37], and the Ising 



model is exactly the two-state Potts model. At p > 3, 
the Potts model shows features, very different from those 
at p = 1 and p — 2. In standard mean-field theory, the 
p > 3 -Potts model has a first-order phase transition in 
contrast to percolation and the Ising model, where the 
phase transitions are of second order. 

Thus, the Potts model provides an essentially more 
wide range of behaviours than the percolation problem 
and the Ising model. The Potts model is related to a 
number of outstanding problems in statistical and math- 
ematical physics, and in graph theory, e.g., the colouring 
graph problem, etc. (for numerous applications of the 
Potts model see Refs. [35,36]). In short, this a basic 
model of statistical mechanics, the direct generalization 
of the percolation problem and the Ising model. Con- 
sequently, before more complicated cooperative models 
on networks, one has to solve the Potts model on the 
simplest complex networks. 

In this paper we present in a detailed form our results 
which we have announced in Ref. [22]. We demonstrate 
a strong effect of a fat tail in the degree distribution of 
a network on a first-order phase transition which occurs 
in the Potts model with p 3. We observe the suppres- 
sion of the first-order phase transition in networks with 
a fat-tailed degree distribution. 

It is convenient to use the power-law degree distribu- 
tion P(k) oc k 1 for parametrization. Then, (k 4 ) di- 
verges for 7 < 5, (k 3 ) diverges for 7 < 4, and (k 2 ) di- 
verges for 7 < 3. 

We find that if p > 3, the first order phase transition 
occurs only for 7 > 3, while for 7 < 3 the phase tran- 
sition is continuous. At 7 = 3 , the phase transition is 
of infinite order and similar to the transition in the Ising 
model on these networks [11]. We also obtain the exact 
solution of the one-state Potts model on networks and 
show that it agrees with the solution of the percolation 
problem on networks [20]. 

In this paper we consider the large networks that have 
a tree- like 'local structure'. In other words, any finite 
environment of a vertex in the infinite network looks like 
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a tree. Interactions are transmitted through edges, from 
vertex to vertex. Thus, it is the distribution of the num- 
ber of connections of a nearest neighbour of a vertex that 
is crucial. In the uncorrelated random networks, this dis- 
tribution is kP(k)/(k). Then the nearest neighbours of a 
vertex have the average number of connections (the aver- 
age degree) {k 2 )/(k), its second nearest neighbours have 
the same average degree, and so on. This value is greater 
than the average number of connections for the entire 
network (k), and it is much greater than (k) if (fc 2 ) is 
large. It is not the mean degree of a network that deter- 
mines its cooperative behaviour but, rather, the average 
degree of the nearest neighbour of a vertex [11]. This 
enables us to estimate a certain characteristic tempera- 
ture of the p-state Potts model on the network using the 
formula T c /J = 1/ ln[(q +p-2)/(q- 2)] for the p-state 
Potts model on a regular Cayley tree with the coordina- 
tion number q = (k 2 ) / (k) [36]. The result is 

J_ (k 2 ) + (p-2) (fc) 

T c (fc 2 ) - 2(fc) ' 1 ! 

This naive estimate is exact for the Ising model [11]. It 
also gives the exact percolation threshold as we will show 
below. Note, however, that the meaning of the critical 
temperature T c is different for p = 1,2 and p ^ 3. At 
p = 1 the parameter q c = cxp(J/T c ) — 1 determines the 
percolation threshold on networks (q is the probability 
that a vertex is retained). At p — 2, T c is the critical 
temperature of the continuous phase transitions in the 
Ising model on networks (with J replaced by 2J). We 
shall show that for p ^ 3 and 7 > 3, the p-state Potts 
model undergoes a first order phase transition, and T c 
given by Eq. 1 is the low-temperature boundary of the 
region of hysteresis phenomena. We shall show, however, 
that if 7 ^ 3 and p ^ 1, the transition in the Potts 
model becomes continuous, and T c is again the critical 
temperature of the transition (see discussion below). 

II. GENERAL APPROACH 

Consider the p-state Potts model with the Hamilto- 
nian: 

n = -jY,s ai , aj -Hj2s ai ,u (2) 

(ij) » 

where the first sum is over all edges of the graph, the 
second one is over all vertices. 5 a ,p = 0, 1 if a 7^ 
and a = 0, respectively. Each vertex i can be in any 
of p states, i.e., a, = 1,2,..., p. We assume a 'ferro- 
magnetic' interaction between the nearest-neighbouring 
vertices, i.e., J > 0. The 'magnetic field' H > distin- 
guishes the state a = 1. Hereafter, we set J = 1. The 
p-state Potts model on the regular Cayley tree is solved 
exactly by using recurrence relations [36]. As networks 
under discussion have a local tree-like structure, we ap- 
ply the method of recurrence relations to our problem. 



Actually, we use the same approach as in our solution of 
the Ising model on networks [11]. 

Consider a vertex with fco adjacent vertices with 
'spins' in states ctij, i = 1, 2, . . . , fc . Due to the local 
tree-like structure, this vertex may be treated as a root 
of a tree. We introduce 

9i,i( a o) = 

+ H^2s ai>1 )/T\. (3) 

{«(} (2m) I 

The labels I and m run only over vertices that belong to 
sub-trees with the root vertex a\^, including this vertex. 
Then the partition function is 

Z^e^n^W- (4) 

Let 

xi,i(a) = gi,i(a)/g lti (l) , (5) 
then the 'magnetic moment' M of the vertex is 

M=^-(s ao , 1 --\ 
p-l\ p/ T 

= 1 (p-l)e g/T -E a#1 II-^M(a) 
" p-l e^ + ^Ylti^M) ' 

where (■■■) T is the thermodynamic average. The param- 
eters x\ t i describe the effects of the nearest neighbours 
on the vertex 0. In turn, are expressed in terms 
of parameters x 2 ,i(a) = g 2 ,i(a)/g2,i{l), 1 = 1,2,..., fc M , 
which describe effects of vertices in the second shell on 
vertices in the first shell, and so on. The following recur- 
rence relation between x n j(a) and x n+ \ t i{a) holds at n 
^ 1 and a ^ 2: 

k — 1 k ■ — 1 

[e H/T + e l/T |j Xn+ll{a)+ J2 "ft *n + M\l 
1=1 P^l,a 1=1 

[e^^ + ^flzn+iAP)}- (7) 

/Mi *=i 

If a vertex n + 1, 1 is a dead end, then x n+ i t i{a) = 1 at 
all a. Deriving the recurrence relations, we started from 
some vertex and then made the recurrence steps along 
sub-trees. While solving the recurrence relations (7), we 
start from distant vertices, i.e., from large n, and descend 
along sub-trees to the vertex. Note that all states with 
the index a ^ 2 are equivalent in respect of statistics. 
Only the state a = 1 is distinguished by the applied field 
H. 
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The recurrence steps Eq. (7) converge exponentially 
quickly to the fixed point which does not depend on a. 



This enables us to put x n j (a) 



in Eq. (7) from the 



very beginning for the sake of simplicity. Then Eqs. (6) 
and (7) take a form: 



M = 



oH/T _ 



n 



feo 



X\ 



e*/ r + (p-l)n£i*i, 



^n,j — y 



n 

i=i 



where we introduce 



y{x) = 



X n +l,l) 



+ {e 1 ' 1 + p-2)x 

e (l+H)/T + ( p _ 1)x 



(8) 



(9) 



(10) 



Note that at H = 0, in the paramagnetic phase x n ,i = 1 
while at H > 0, x n ,i ^ 1- We stress that Eqs. (8) and 
(9) are exact for any tree-like graph. 



III. EXACT SOLUTION AT P = 1: 
PERCOLATION ON NETWORKS 

The one-state Potts model on networks at H = is 
of a special interest, as it relates to percolation on net- 
works. This limiting case can be solved exactly for an 
uncorrelated random graph with an arbitrary degree dis- 
tribution. For this let us consider the recurrence relation 
(9) at p = 1 and H = 0: 



IV. FIRST ORDER PHASE TRANSITION 



It is well known that the p-state Potts model, in the 
framework of the standard mean-field theory undergoes 
a first order phase transition for all p ^ 3 [35,38]. The 
general approach derived above enables us to consider, 
at H = 0, the influence of the fat tail of the degree dis- 
tribution on the transition. 

At p ^ 3, in order to average over the ensemble of ran- 
dom graphs with a given degree distribution we use the 
effective medium approach developed in Ref. [11] for the 
Ising model on networks. First at all, as at H ^ we 
have x n j ^ 0, it is convenient to introduce positive pa- 
rameters h n ,i- x n j = cxp(—h n j). They are independent 
random parameters and may be considered as random 
effective fields acting on a vertex in the n-th shell from 
neighbouring 'spins' in the n + 1-th shell. Then, Eqs. (8) 
and (9) take the form: 



M = 



s g/T -exp(-E£ifei,0 
e"/T + (p-l)exp(-££i>*i,0 



K,j = - In < y 



exp 



+i,i 



i=i 



(14) 
(15) 



e- 1 / T + (l-e- 1 ' T ) 



i n +i,i 



(11) 



i=i 



Note that in an uncorrelated random graph, the degrees 
k n j are independent random variables. Since x n+ u in 
this equation depends only on k m j with m n + 1, one 
can average the left and right sides of the equation over 
the ensemble of random graphs with a degree distribution 
P{k) and find the self-consistent equation for the average 
(x n j). In the limit n — > oo, the fixed point 
the exact recurrence relation is given by 

q ^P(fc)fca 



x = 1 - q + 



„fe-i 



(k) 



) — > x of 



(12) 



where we introduce the parameter 



9 = 1 



,-l/T 



(13) 



It is important to note that Eq. (13) establishes the re- 
lation between the one-state Potts model and site perco- 
lation on uncorrelated networks. The latter is described 
by Eq. (12) [20]. The parameter q has the meaning of 
the retained fraction of vertices. The critical tempera- 
ture T c , Eq. (1), at p = 1, determines the percolation 
threshold q c = (k) /((fc 2 ) — (k)). The strong influence of 
the fat tail of the degree distribution on percolation crit- 
ical exponents has been revealed and studied in detail in 
Ref. [20]. 



At dead ends we have = 0. At H = in the 

paramagnetic phase h n j — 0, while in the ordered phase 
h n ,i 7^ 0. Equations (14) and (15) determine the magne- 
tization M of a graph as a function of T and H. 

While solving the recurrence relations (15), we start 
from distant spins with h w and descend along sub- 
trees to a. In the limit n — > oo, the parameter is the 
fixed point of the recurrence steps. The thermodynamic 
behaviour is determined by this fixed point. 

The right-hand sides of Eqs. (14) and (15) depend 
only on the sum of the independent and equivalent ran- 
dom variables h n j. So, let us use the following ansatz 
[11]: 



Y.h^i-kh + Oik 1 / 2 ), 



(16) 



where h = (h n i) is the average value of the 'effective 
field' acting on a vertex. This approximation takes into 
account the most 'dangerous' highly connected vertices 
in the best way. With this ansatz, averaging over the 
ensemble of random graphs and applying the ansatz (16) 
to Eqs. (14) and (15), we obtain 
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(M) =J2 P (k)- 



oH/T 



-kh 



(17) 



e H / T + (p - l)e- kh ' 

h = - {ky 1 ^P{k)k\ny[e- {k - 1 ^ h ]= G{h) . (18) 
fe 

The parameter ft plays the role of the order parameter. 
At H = 0, ft = above T c and ft, > below T c . 

Let us describe the thermodynamic properties of the 
Potts model with p ^ 3 on the infinite networks at H = 0. 
For this, one must solve the equation of state (18). 

At first, we consider the case (k A ) < oo. The character 
of the transition and a characteristic temperature may 
be found from the analysis of the expansion of G(h) over 
small ft. If p ^ 3 it is enough to take into account only 
first two terms: G{h) = g\h + g 2 h 2 + ... where 



that is expected for the infinite order phase transition at 
7 = 3 [11]. 



5i = 



92 = 



(k(k-l))(e 



l/T 



1) 



(jfe)(eVr +p _i) ' 

(fc(fc-l) 2 )(e 1/r -l)(p- 
2(fc) (eVT+p- 1)2 



2) 



(19) 
(20) 



One sees that at high temperatures, the coefficient g\ w 
0. When temperature decreases, g\ increases and be- 
comes larger than 1. The point where gi = 1 is the spe- 
cific point. The corresponding temperature T c is given by 
Eq. (1). As <72 > at all temperatures, the Potts model 
with p ^ 3 undergoes a first order phase transition from 
the paramagnetic phase to the ordered one. Note that 
for the Ising model (p = 2) the coefficients g 2 = and 
gs < 0, and the phase transition is continuous [11]. In the 
paramagnetic phase the probabilities to find any vertex 
in the states a = 1,2, . . . ,p are equal, which corresponds 
to the solution ft = 0. In the ordered phase, one of the 
states, a = 1 in our consideration, has a larger proba- 
bility in comparison to other states with a ^ 2. This 
corresponds to the solution with ft > 0. 

In order to find the temperature behaviour of the mag- 
netization in the first order phase transition, we have 
solved numerically the Eq. (18). At large 7, the system 
undergoes a standard first order phase transition (see 
the temperature behaviour of M in Fig. 1 for 7 = 6). 
The decrease of 7 modifies the second term in the ex- 
pansion of G(h): at 7 = 4 we have (k A ) — > 00 and 
G(ft) = g x h + g' 2 h 2 ln(l/ft) + . . ., at 3 < 7 < 4 wc find 
G{h) = gih + g' 2 h'^ 2 + ... In the range 3 < 7 < 4, the 
derivative dG/dh at small ft remains to be positive and 
the transition is still of the first order (sec in Fig. 1 our 
results for 7 = 4 and 3.2). Results of the numerical solu- 
tion of Eqs. (18) and (17) reveals that when 7 approaches 
3 from above, the jump of the magnetic moment in the 
first order phase transition tends to zero. Moreover, it is 
interesting to note that the temperature behaviour of M 
in a wide temperature range just after the jump follows 
the exponential law M ~ exp(— cT) with the constant c 
which depends on the complete degree distribution P(k) 
(see the results in Fig. 1 for 7 = 3.2). It is the behaviour 




FIG. 1. The temperature dependence of the magnetization 
of the p-state Potts model with p = 5 on uncorrelated ran- 
dom networks with the degree distribution P(k) = Ak~"' for 
different exponents 7: 1) 7 = 6; 2) 7 = 4; 3) 7 = 3.2. The 
temperature region between dotted lines for each curve is the 
hysteresis region. The up and down arrows correspond cool- 
ing and heating, respectively. 

In the case of the first order phase transition, there is 
a temperature region, T2 < T < Ti, of hysteresis phe- 
nomena owing to the existence of metastable states. For 
the Potts model, the low temperature boundary T2 is the 
temperature below which the ordered state with ft ^ 
is the only stable state. In other words, below T 2 the 
free energy as a function of the order parameter has only 
one minimum corresponding to ft 7^ 0, while the solution 
ft = corresponds to the maximum. One can show that 
T 2 is determined by the equation g\ — 1, i.e. T 2 = T c . In 
the range T 2 < T < T lt the states with ft ^ and ft = 
correspond to the minimum of the free energy, i.e. one of 
these states, having larger energy, is metastable. At tem- 
peratures T > Ti the paramagnetic state, ft = 0, is the 
only possible solution of the equation of state (18). T\ 
is determined by the set of Eqs. (18) and dG(h)/dh = 1. 
The hysteresis region is shown in Fig. 1 at different 7. 
One can see that its width AT = T\ — T 2 increases with 
increasing 7, however, AT/T 2 — > 0. 



V. CONTINUOUS TRANSITION FOR 2< 7 ^3 

The case 7 = 3. — Here, the second moment (& 2 ) 
diverges. Using, for brevity, the continuum approxi- 
mation for the degree distribution, we obtain G(h) w 
((k) h/(pT))\n[p/({k)h)}. One sees that at small ft, the 
second derivative d 2 G(h)/dh 2 is negative in contrast to 
the case 7 > 3 where d 2 G(h)/dh 2 is positive. It means 
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the change of the order of the phase transition. Instead 
of the first order phase transition discussed above, the p- 
state Potts model with p ^ 3 undergoes an infinite order 
phase transition at the critical temperature T Cl Eq. (1), 
similarly to the Ising model [11] and percolation [20]. 
This conclusion also agrees with the phcnomcnological 
theory of critical phenomena in complex networks [22]. 

One should emphasize that when (fc 2 ) diverges, the 
critical temperature T c is infinite for the infinite net- 
works [see Eq. (1)]. However, in any finite network, 
(fc 2 ) < oo, and T c is finite, although it may be very high, 
T c = (k 2 )/((k)p) (see below). At temperatures, which 
are much less than T c , but where ft < 1, so T > 1, we 
obtain 

h = (p/(k))e- 2 'P T /^ , M £* e - 2 P T / < fc > , (21) 

Without the continuum approximation, we have, instead 
of (k) in the exponential, a constant which is determined 
by the complete form of P(k). 

The case 2 < 7 < 3. — Again T c for large networks is 
very high and the phase transition is continuous. Using 
the expansion G(h) = g((k) /(pT))h,i- 2 , we find, in the 
range 1 <C T <C T c , that the Potts model demonstrates 
the behaviour 

h,M ~ 7-1/(8-7) (22) 

which is quite similar to the Ising model [11]. 

In accordance with Eq. (1), T c diverges when (fc 2 ) — > 
00. However, in finite networks, (fc 2 ) is finite because 
of the finite-size cutoff of the degree distribution. In 
scale-free networks, it is usually estimated as k cu t ~ 
fcoA 1 ^ 7-1 ), where N is the total number of vertices in a 
network, fco is a 'minimal degree' or the lower boundary 
of the power-law dependence, and (fc) « fco(7 — 1)/(7 — 2). 
Then, using estimates from Refs. [4,20,29] we obtain 

T c « 1-i at 7 = 3 , 

V 

T c - J (7 7, 2)2 - (k)N( 3 -^/^-V for2< 7 <3. (23) 
p(3-7)(7-l) 

These expressions generalize the finite-size effect ob- 
tained for the Ising model on networks [11]. 

VI. DISCUSSION 

Our investigations of the Potts and Ising models on 
uncorrelated random networks demonstrate a general 
strong effect of fat tails in the degree distribution on the 
order of the phase transition and its critical tempera- 
ture. In the Potts model with p ^ 3 states, the phase 
transition of the first order occurs only when the second 
moment (fc 2 ) is finite, i.e., at 7 > 3. When (fc 2 ) diverges 
the behaviour of the Potts and Ising models arc similar. 
They undergoes the infinite order phase transition. Wc 



suggest that this phenomenon (the suppression of a first 
order phase transition in favour of an infinite order one) 
is of a general nature and it takes place in other coop- 
erative models with short range interaction on random 
networks. 

The phase diagram of the Potts model on uncorrelated 
networks was also studied in Ref. [23] in the framework 
of a simple mean field approach. One should note that 
our results essentially differ from those in Ref. [23]. The 
reason for this difference is evidently the simplified mean 
field theory of Ref. [23]. 
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